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Pion scattering on the lattice with chirally improved fermions* 
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We report preliminary results for the 7 = 2 pion scattering length ao, calculated with chirally improved fermions. 
After chiral extrapolation our results are in good agreement with both, experimental and theoretical predictions. 
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1. Introduction 

Since the pion is the goldstone boson of the 
spontaneous breaking of chiral symmetry, observ- 
ables related to the pion give insight into the chi- 
ral dynamics of QCD. For example the pion scat- 
tering length is a parameter entering the effective 
chiral lagrangian of QCD and a precise determi- 
nation of this observable in an ab-initio lattice 
calculation is desirable. On the lattice scattering 
data of stable particles can be calculated using 
Liischer's formula 1 which relates the scattering 
length to the energy W of the scattering state. 
The crucial technical challenge is a very precise 
determination of W and advanced strategies for 
error reduction have to be used. For chirally 
sensitive quantities, such as the pion scattering 
length, it is important to be able to work with 
small pion masses. Using fermions based on the 
Ginsparg- Wilson equation it is now possible to 
reach pion masses much lower than what was ac- 
cessible with more traditional formulations. This 
implies that the chiral extrapolation of the results 
to the physical pion mass becomes more reliable. 

2. Setting of our calculation 

Our calculation is based on the chirally im- 
proved (CI) Dirac operator CI fermions are 
a systematic approximation of a solution of the 
Ginsparg- Wilson equation and have been tested 
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conf. 
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7.90 


0.15 


2.4 


16^x32 


97 


B 


7.90 


0.15 


1.8 


12^x24 


100 


C 


8.35 


0.10 


1.6 


16^x32 


91 



Table 1. Parameters and statistics of our quench- 
ed ensembles. 

successfully |3] in spectroscopy down to pion 
masses of m^r = 265 MeV. 

We use gauge configurations generated with the 
Liischer-Weisz action. Our quark masses are in 
the range of am = 0.015 to am = 0.2. The scale 
was determined in 0] using the Sommer param- 
eter. The parameters of the three ensembles we 
used are listed in Table 1. 

3. Liischer's formula 

If two stable particles are confined in a box of 
size L, the energy W of the corresponding scat- 
tering state can be expanded |1| in a power series 
in L^^. The coefficients of this series are related 
to the elastic scattering amplitude, 
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* Poster presented at LATTICE 2004. This work was sup- 
ported by DFG and BMBF. 



where the omitted terms are 0{L ^) and ci, C2 
are numerical constants that are known analyti- 
cally (ci = -2.837297, ca = 6.375183). Thus, for 
calculating the s-wave scattering length oq, one 
has to compute and W on the lattice and to 
solve the cubic equation ^ for ao . The challenge 
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is to obtain the difference W — 2-1X1^^ with high ac- 
curacy since the solution of the cubic equation is 
very sensitive to this difference. 

4. Two pion scattering state energy 

We calculate the pion mass m-^ from a 2- 
parameter fit of the pseudoscalar correlator (T 
denotes the temporal extent of the lattice) 

= (P(t)pt(o)} ^ Acosh(TO,[r/2-t]),(2) 

where generates I — 1 pseudoscalar mesons 
and P denotes the zero-momentum Fourier trans- 
form of the corresponding annihilator. 

The energy W of the scattering states can be 
extracted from 



{C,{t)) = (P(t)P(t)pt(o)pt(o)) 



(3) 



In both, the 2- and the 4-point correlators, we use 
Jacobi smeared quark sources and sinks. We also 
experimented with wall-type sources but found 
no improvement of the signal. 

While the pion mass can be extracted quite ac- 
curately from |5J, a direct determination of the 
scattering state energy W from a fit to © is un- 
satisfactory. A successful method 5 to reduce 
the statistical noise is to consider the ratio 
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where fluctuations of the correlators on single 
configurations cancel. We first determine 
from and subsequently use this value in the 
3-parameter fit l@J which determines W . Our fits 
are fully correlated and errors were determined 
with the jackknife method. 

5. Chiral extrapolation 

For the chiral extrapolation of our data we use 
chiral perturbation theory. In one-loop chiral per- 
turbation theory one finds '5 
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Ga is a (scale dependent) combination of LEC's 
and /o is the pion decay constant in the chiral 
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Figure 1. Results for Pi(mJ) = aa/mT^ as a func- 
tion of the squared pion mass. The open symbols 
are our data, the filled symbols their extrapola- 
tion to the physical pion mass. The curves repre- 
sent our fits. The experimental value is indicated 
by the asterisk. 



limit. However, we do not fit our data directly to 
the form , but consider two combinations of oq 
with m.^.. To remove the leading m7r-dependence 
of the scattering length we study the ratio 
ao 



F,{ml) 

We also consider the product 
F2{ml) = aom^ , 



(6) 



(7) 

which is often used by experiments. We perform 
the chiral extrapolation for both Fi and F2 using 
the accordingly modified right hand side of 
We furthermore experimented with adding higher 
order terms from chiral perturbation theory but 
did not observe a large effect. 

6. Results 

In Fig. 1 we present our results for the ratio Fi. 
The open symbols are our data and the filled sym- 
bols arc their extrapolation to the physical pion 



3 



° 16^32, L=2.4 fm (ens. A) 

* 12^24, L=1.8fm(ens. B) - 

o 16^32, L=1.6fm (ens. O- 

[^;v ■)]<- experiment 




-2.0 - 



_2 5LI , 1 , 1 , 1 , 1 

■ 0.0 0.5 1.0 1.5 2.0 

ml [GeV'] 

Figure 2. Same as Figure 1, but now for the prod- 
uct F2(m^) — ao m^r. 



mass. The asterisk represents the experimental 
value. When comparing the two ensembles B (tri- 
angles) and C (circles) which are similar in size 
(L = 1.8 fm and L = 1.6 fm) but have a different 
lattice spacing (a — 0.15 fm and a — 0.1 fm) we 
find, at least for larger pion masses, good agree- 
ment of the data points. This indicates that dis- 
cretization errors are small, similar to what was 
already found in quenched hadron spectroscopy 
with the chirally improved Dirac operator 0. A 
comparison with the data obtained on our larger 
volume (ensemble A, L — 2.4 fm, a = 0.15 fm) 
shows, however, that at the smaller volume with 
L = 1.8 fm we still observe a finite size effect. 

The same picture also holds for the analysis of 
F2. Again the two ensembles with small physical 
volume but different lattice spacing agree reason- 
ably well, indicating that scaling violations are 
small. The larger volume gives lower values for 
F2, showing that at L = 1.8 fm finite size effects 
cannot be neglected. 

Because of the good scaling and the finite size 
effect we expect the results from ensemble A 
(16^ X 32, a = 0.15 fm, L = 2.4 fm) to be our best 



extrap. our data experiment X"PT 
Fi(ml) -0.0453(35) -0.0454(34) -0.0444(10) 
F2iml) -0.0425(97) -0.0454(34) -0.0444(10) 

Table 2. Our final results for aom^ compared to 
experimental data and to 2-loop chiral perturba- 
tion theory. 

estimate. For this ensemble we list in Table 2 our 
results for ooto^ obtained from the extrapolation 
of both Fi(toJ) and ^2(771^) to the physical pion 
mass. We compare our final results to recent ex- 
perimental data 7; and to a two-loop result from 
chiral perturbation theory |S| . Given the fact that 
our calculation is quenched, the agreement of our 
results with the experimental number is surpris- 
ingly good. 
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